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Abstract. The filamentary structure of the early universe has until now only been seen in numerical simulations. 
Despite this lack of direct observational evidence, the prediction of early filamentary structure formation in a 
Cold Dark Matter dominated universe has become a paradigm for our understanding of galaxy assembly at high 
redshifts. Clearly observational confirmation is required. Lyman Break galaxies are too rare to be used as tracers 
of filaments and we argue that to map out filaments in the high z universe, one will need to identify classes of 
objects fainter than those currently accessible via the Lyman Break technique. Objects selected via their Lya 
emission, and/or as DLA absorbers, populate the faintest accessible part of the high redshift galaxy luminosity 
function, and as such make up good candidates for objects which will map out high redshift filaments. Here we 
present the first direct detection of a filament (at z=3.04) mapped by those classes of objects. The observations are 
the deepest yet to have been done in Lya imaging at high redshift, and they reveal a single string of proto-galaxies 
spanning about 5 Mpc (20 Mpc comoving). Expanding the cosmological test proposed by Alcock & Paczyhski 
(1979), we outline how observations of this type can be used to determine at z=3. 
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1. Introduction 

■ For the past three decades, following the first prediction of 
pancakes by Zel'dovich 1970, simulations have been ahead 
of the observations when it comes to describing the first 
structures to form at high redshift. Numerical simulations 
of structure formation based on the Cold Dark Matter sce- 
nario predict that the first large scale structures to form 
are voids and filaments (Klypin & Shandarin 1983; White 
et al. 1987; Evrard et al. 1994; Rauch et al. 1997). In this 
picture any sightline through the early universe will inter- 
sect a large number of voids and filaments, and the varia- 
tion in the density of neutral hydrogen along the sightline 
is then observed as the Lyman forest (Petitjean et al. 1995; 
Theuns et al. 1998). Growth of density fluctuations along 
the filaments will lead to formation of lumps of cold, self- 
shielding gas and those are identified, in the simulations, 
as regions of star formation (Katz et al. 1996; Haehnelt et 
al. 2000). Because of the high column density of neutral 
hydrogen in such clouds, they are identified observation- 
ally as strong absorbers and the strongest as Damped Lya 
Absorbers (DLAs). By poking random sightlines through 
a virtual universe one may simulate observations, and a 
given model universe will hence predict a specific correla- 
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tion between DLA systems and the galaxies hosting the 
DLAs (Katz et al. 1996). Comparison to real observations 
of DLA galaxies (M0ller & Warren 1998) has shown that 
there is very good agreement between observations and 
simulations. This agreement is encouraging, but it would 
be of great interest if one could observationally map out 
the actual high redshift filaments. Apart from being a fun- 
damental test of the validity of the simulations and the 
structure formation paradigm they represent, knowing the 
scale size of the filaments at different redshifts will help 
constrain both cosmological parameters and the mecha- 
nisms of galaxy formation. 

Unfortunately such a map cannot be constructed via 
absorption studies, because there is no sufficiently tight 
mesh of background quasars available. The only way to 
proceed is to attempt to find enough centers of star for- 
mation to map out the filament in its own light, as has 
been done at low redshift (z<0.2) using galaxy surveys 
(Gregory & Thompson 1978; de Lapparent et al. 1991; 
Bharadwaj et al. 2000). At low redshifts filaments remain, 
but many of the galaxies have drained into the filament 
nodes where todays galaxy clusters are formed. In order to 
identify objects in a high redshift filament one might first 
try a search for Lyman Break Galaxies (LBGs, Steidel 
& Hamilton, 1992). Unfortunately only relatively bright 
galaxies can be found with this technique, and sparse sam- 
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Table 1. Redshifts and positions of seven Lya emitters 
and a Lya absorber in the field of Q1205-30. The positions 
are given relative to the quasar coordinates: 12:08:12.7, - 
30:31:06.10 (J2000.0). The uncertainty on the redshifts is 
0.0012 {la). 



N 



Object 


ARA (arcsec) 


Adecl. (arcsec) 


redshift 


S7 


-143.3±0.6 


41.9±0.2 


3.0402 


S8 


-141.5±0.6 


59.7±0.2 


3.0398 


S9 


-124.6±0.5 


63.4±0.2 


3.0350 


S10 


-119.9±0.5 


59.8±0.2 


3.0353 


Sll 


-77.8±0.3 


0.9±0.1 


3.0312 


S12 


-43.9±0.2 


54.4±0.2 


3.0333 


S13 


68.3±0.3 


-52.1±0.2 


3.0228 


abs 


0.0 


0.0 


3.0322 



pling of the filamentary structure does not allow it to be 
seen easily. However, both DLA galaxies and galaxies se- 
lected for their Lya emission, sample the high redshift 
galaxy population at much fainter magnitudes than do 
the LBGs, and should therefore give a better sampling of 
the high redshift structure (Fynbo et al. 1999; Haehnelt et 
al. 2000; Fynbo et al. 2001 (FMT2001)). This has recently 
been independently confirmed, by deep narrow band Lya 
imaging in a field of known LBGs, which reveals about 
a factor of 6 more candidate Lya galaxies than LBGs 
(Stcidcl et al. 2000). A hint that filamentary structure is 
indeed present, comes from a few cases where narrow-band 
Lya imaging of DLAs was attempted, and one or two ad- 
ditional neighbour Lya emitters were found. In particular 
it has been demonstrated that when three DLA and/or 
Lya sources are detected in a field, they show a marked 
tendency to be aligned (M0ller & Warren 1998). 

2. Objects in the field of Q1205-30 

In Fynbo et al. 2000 we reported on a new very deep, 
1.1 x 10~ 17 erg s -1 cm ~ 2 (5er), narrow band search for 
Lya emitters in a 28 arcmin 2 field centred on the quasar 
Q 1205 30. The filter was tuned to the wavelength of a 
strong Lya absorption line in the quasar spectrum. On 
March 4-5, 2000 we used the ESO Very Large Telescope 
on Paranal to obtain spectra of the candidate Lya emit- 
ters in the field. Details of the spectroscopic observations 
and data reductions have been submitted for publication 
(FMT2001), but the results are summarised in Table 1 
where we list the position and redshift of each of the seven 
confirmed Lya emitting objects. In the table we also give 
the coordinates and redshift of the absorber. 

In Fig. 1 we show a sky-coordinate map of the field. 
While there may here be a hint of an elongated structure 
running from the SE corner of the CCD frame towards 
the NW, the statistical evidence for a non-random dis- 
tribution is weak. However, this may be because we are 
looking almost along a filament. If this was the case, we 
would expect the redshifts to change monotonically along 
the structure suggested by the dotted line in Fig. 1. This 
should give rise to a strong correlation between redshift 
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Fig. 1. Position of the eight objects on the CCD. In 
this and all following figures solid hexagons mark the Lya 
emitters, the open hexagon marks the absorber. 
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Fig. 2. Left: Redshift difference plotted against relative 
distance on the sky (in CCD pixels). The correlation coef- 
ficient is r s = 0.73. In random distributions the correlation 
parameter is worse than that in 99.7% of all cases. Right: 
Redshift plotted against projected ID distance (positions 
projected onto the dotted line in Fig. 1). 



and position on the CCD. We performed two tests to as- 
sess this correlation, one Monte Carlo and one analytical: 
1) For each of the eight objects listed in Table 1 we 
measured the projected distance on the sky, as well as 
the difference in redshift, to each of the other seven ob- 
jects. Those measurements are reproduced in Fig. 2 (left). 
A Spearman correlation test of the data points gives a 
correlation coefficient of r s = 0.73. We then produced 
3 x 10 6 realizations of a random distribution of eight ob- 
jects in the observed volume, keeping the position and 
redshift of the absorber fixed. In only 0.3% of the simula- 
tions did the same test give a similar or better correlation. 
A more fair comparison is probably to compare against a 
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Fig. 3. 3D distribution of the eight objects seen from 4 different viewing-anglcs. fn each of the figures the 3D arrow 
points in our viewing direction into the sky, and the spiral pattern maps out a cylinder with diameter 400 kpc (see 
Table 2). The box marks the volume of space observed with our narrow-band Lya filter. Top left: Here we have rotated 
the view to look along the filament. The thick lines mark the front "entrance window" of the box (corresponding to 
our CCD image). Top right: The box is here rotated 90 degrees to the right, hence viewing the filament from the 
left side compared to the end-on view. Bottom right: Same as top right but rotated 90 degrees around the filament 
to give a view of the filament as seen from "above" the view in top left. Bottom left: View from a random angle to 
give an impression of the 3D structure. 



random distribution with an imposed two-point correla- 
tion function as reported by Giavalisco & Dickinson 200 f 
(£(r) = (r/r )~ 7 , r = Ll/i^Mpc, 7 = 2.2). In this case 
0.4% of the simulations gave a similar or better correla- 
tion. 

2) Encouraged by this tentative suggestion of a fila- 
ment, we proceeded to project the position of each object 
on the CCD onto the line drawn in Fig. 1. In Fig. 2 (right) 
we plot the redshift of each object against the projected 
position. A Spearman correlation test now gives r s = 
0.881 for eight objects. In this case we are dealing with 
statistically independent data points, which confirms the 
correlation with a probability of 99.82%. 

These tests confirm that the eight objects form a tight 
string-like structure in redshift space. Assuming that the 
rcdshifts are caused by Hubble flow, this structure takes 
the form of a 3D filament similar to those in numerical sim- 



ulations. The objects are fainter and much more common 
than the typical LBGs and they presumably represent an 
earlier evolutionary state than the brighter LBG galaxies. 
Evidence for this is found in the fact that only 20% of the 
LBG population has Lya in emission, and as a class they 
show evidence for internal dust absorption suggesting high 
metallicity (Steidel et al. 2000). In contrast, the numerous 
faint objects with Lya in emission have much lower star 
formation rates than the LBGs and are essentially dust 
free (FMT2001). In the simulations the final fate of sev- 
eral such objects is to merge into one single galaxy, and for 
this reason we shall here refer to them as "Lya Emitting 
Galaxy-building Objects" (LEGOs). 

3. The Q1205-30 filament at z=3.04 

In Fig. 3 we show the 3D structure of the volume of space 
we are sampling from 4 different viewing- angles. Our sam- 
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Table 2. Filament parameters for different cosmologies. 
In all cases have we used H Q = 65 km s _1 Mpc -1 . We 
denned the length as the distance between the two out- 
hermost LEGOs. 





Model 




rms (kpc) 


Radius (kpc) 


Length (kpc) 


Q 


= 0.2, Q A 


= 0.0 


280 


400 


4400 


Q 


= t.o, n A 


= 0.0 


180 


260 


2800 


Q 


= 0.3, £l A 


= 0.7 


280 


400 


4800 


Q 


= 0.1, Qa 


= 0.9 


390 


560 


7700 



pie inside this volume is complete down to a limiting Lya 
flux of 1.1 x 1CT 17 erg s _1 cm ~ 2 (5<r). The redshift dis- 
tribution of the LEGOs is consistent with a random dis- 
tribution when one allows for the selection defined by the 
filter transmission curve (FMT2001), and we hence con- 
clude that the distribution along the filament is consistent 
with a random distribution. In other words, we have not 
detected the ends of the filament, and it may extend fur- 
ther in either direction. 

In each of the four figures in Fig. 3 we have enclosed the 
LEGOs in a hollow "tube" which is meant to illustrate the 
surface of the filament. The same filament-surface is shown 
in all the figures. In particular in Fig. 3(top left) where the 
filament is viewed from the end, the tube takes the form of 
a circle. To define the scales of the filament we have chosen 
three parameters: The length, the radius of the tube and 
the rms of the projected distribution of the objects around 
the center of the filament when the filament is viewed from 
the end. The actual values of those parameters will depend 
on the chosen cosmology. In Table 2 we list the results 
for a representative selection of cosmologies. Note that to 
represent the length we have chosen the largest distance 
between two detected LEGOs. As discussed above this 
is in reality a lower limit to the length as we have not 
yet determined where the endpoints are. The comoving 
length is hence found to be minimum 20 Mpc, in good 
agreement with recent predictions (20-25/i , Dcmiahski 
& Doroshkevich 1999). 

In an O = 0.2, Ha = 0.0 cosmology the volume sam- 
pled by our observations is 30 Mpc 3 . The filament fills 
only about 7% of this volume, consistent with a picture of 
the early universe where the non-linear processes leading 
to the formation of proto-galaxies are confined within thin 
filaments. 



4. High z filaments as cosmological probes 

It is interesting to note that five of the seven LEGOs 
remain undetected in the individual broad-band images 
(FMT2001), which reach limiting continuum AB magni- 
tudes of 1=26.9 and B=27.7 (2a upper limits). Our argu- 
ment at the outset that searches for LBGs will not detect 
this type of filament is clearly vindicated. Moreover most 
of the LEGOs reported here have Lya fluxes very close to 
our detection limit. Hence, a slightly shallower Lya search 
would likewise have failed to detect the filament. Finally, 



had the orientation of the filament not been so closely 
aligned with our line of sight, then only a fraction of the 
LEGOs would have fallen within the CCD field of view, 
again making a detection impossible. In this particular 
case only the fortunate near-alignment and the close-to- 
critical detection limit made the detection possible. Lya 
searches to the same or slightly deeper detection limits but 
over larger fields should easily pick out several such fila- 
ments, thereby making it possible to map the filamentary 
structure of a significant volume of the early universe. 

This opens up an interesting perspective as it will pro- 
vide an independent way to determine the cosmological 
constant Oa- The conversion of the observed volume into 
proper coordinates depends on the chosen cosmological 
model. The effect of using a non-zero ft a is to significantly 
stretch the volume along the sightline, thereby causing the 
"length to radius ratio" of filaments to be different for 
filaments seen from the side and filaments seen end-on. 
This is a new realisation of the test proposed by Alcock & 
Paczyfiski 1979, who considered the elongation of an ideal- 
ized spherically symmetric distribution of emitters. In the 
case of emitters distributed in filaments the test becomes 
somewhat more complex because of the lack of spheri- 
cal symmetry, and one will need to consider the statis- 
tics of a sample of filaments. At the same time, however, 
the lack of spherical symmetry offers several independent 
tests. We consider three tests to be of practical use. First 
one may consider the distribution of inclination angles of 
filaments which should be isotropic, secondly the ratio of 
the rms to the length (as given in Table 2) which must be 
independent of filament inclination angle, and third the 
cross-section of the filaments which must in the mean be 
circular for all inclination angles. In practice one should 
determine the value of CIa which provides the best fit to 
all three tests simultaneously. As already pointed out by 
Alcock & Paczynski, random motion of objects with re- 
spect to the Hubblc-flow will add some noise to any such 
test. We shall return to a detailed discussion of this in a 
forthcoming paper. 

This new set of Oa tests is useful in the redshift range 
z = 2-4, is within the capability of current instrumenta- 
tion and is independent of the successful SN (Riess et al. 
1998; Perlmutter et al. 1999) and weak-lensing (Maoli et 
al. 2001) cosmological tests which currently can be used 
out to z « 1. With the complement of large telescopes 
currently available, we predict rapid progress in this field. 
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